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We report on a systematic study of the Coulomb blockade effects in nanofabricated narrow con¬ 
strictions in thin (Ga,Mn)As films. Different low-temperature transport regimes have been observed 
for decreasing constriction sizes: the ohmic, the single electron tunnelling (SET) and a completely 
insulating regime. In the SET, complex stability diagrams with nested Coulomb diamonds and 
anomalous conductance suppression in the vicinity of charge degeneracy points have been observed. 
We rationalize these observations in the SET with a double ferromagnetic island model coupled 
to ferromagnetic leads. Its transport characteristics are analyzed in terms of a modified orthodox 
theory of Coulomb blockade which takes into account the energy dependence of the density of states 
in the metallic islands. 
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I. INTRODUCTION 

(Ga,Mn)As, discovered by Ohno et al. 1 nearly two 
decades ago, is by now the best studied ferromagnetic 
semiconductor®®. An interesting aspect of this mate¬ 
rial are large magnetoresistance effects which were dis¬ 
covered in nanofabricated narrow constrictions in thin 
(Ga,Mn)As filmJ^H® While the effects were initially in¬ 
terpreted in terms of the tunneling magnetoresistance 
(TMRj^l and tunneling anisotropic magnetoresistance 
(TAMR^, it was proven later that the interplay with 
Coulomb bl ockad e is also relevant in narrow (Ga,Mn)As 
constrictions 101 ’®. The origin of this Coulomb blockade 
anisotropic magnetoresistance (CBAMR) effect are sub¬ 
stantial nanoscale fluctuations in the hole densitj® form¬ 
ing puddles of high hole density separated by low con¬ 
ducting regions. (Ga,Mn)As is known to be a strongly 
disordered material. Its hole density is close to the metal- 
insulator-transition. Little variations in the hole density 
caused by local potential fluctuations can lead to an in¬ 
trinsic structure consisting of metallic islands separated 
by insulating areas. It was shown that the magnetoresis¬ 
tance depends, in the presence of Coulomb blockade, not 
only on an applied gate voltage but can also be tun ed by 
changing the direction of the applied magnetic field 10 11 . 
The latter results from the dependence of the Fermi en¬ 
ergy on changes in the magnetization SM and was mod¬ 
eled phenomenologically by Wunderlich et al!®. If trans¬ 
port occurs through a narrow nanoconstriction, single 
electron tunneling (SET) between islands of high carrier 
density becomes relevant. Thus it is not surprising that 
the bias and temperature dependence of the magnetore¬ 
sistance for different magnetization directions could be 
fitted with a model for granular metals in which metallic 
islands are separated by insulating region^ 11 ! Because 
of the nanoscale size of the involved ’’metallic” islands, 
the Coulomb-charging energy U is the dominating energy 
for transport across the nanoconstriction at low tempera¬ 
tures and small bias voltages 14- Since usually more than 
one island is involved in transport, Coulomb blockade di¬ 


amonds, where the resistance is plotted as a function of 
both bias and gate voltage, revealed a very complex and 
irregular pattern. Up to now a detailed experimental 
and theoretical analysis of the Coulomb blockade effects 
in (Ga,Mn)As nanoconstrictions in the single-electron- 
transistor regime is still missing. 

The aim of this work is a systematic study of the 
Coulomb blockade effects in nanofabricated narrow con¬ 
strictions in thin (Ga,Mn)As films. By means of a two 
step electron beam lithography (EBL) technique we fab¬ 
ricated well defined nanoconstrictions (NC) of differ¬ 
ent sizes. Depending on channel width and length, for 
a specific material, different low-temperature transport 
regimes could be observed, namely the ohmic regime, 
the single electron tunnelling regime and a completely 
insulating regime. In the SET regime, complex stability 
diagrams with nested Coulomb diamonds and anoma¬ 
lous conductance suppression in the vicinity of charge 
degeneracy points have been observed. In order to un¬ 
derstand these observations we propose, for a specific 
nanoconstriction, a model consisting of two ferromag¬ 
netic islands coupled to ferromagnetic leads. We study 
its transport characteristics within a modified orthodox 
theory of Coulomb blockade which takes into account the 
energy dependence of the density of states in the metallic 
islands. 

The paper is structured as follows: Sect. [D] explains 
the fabrication process of the samples. In Sect. |III| 
the measurement setup is presented. The next section, 
Sect. |IV| summarizes the results of the measurements, 
giving a first interpretation in terms of a double island 
structure within a classical orthodox model of Coulomb 
blockade 12 17 . In Sect, [v] we present the details of the 
ferromagnetic double island model, study its transport 
characteristics and make a direct comparison with the ex¬ 
perimental results in Sect. |VI A| Conclusions are drawn 
in Sect. IVIII 
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II. SAMPLE FABRICATION 

Our NC-devices were fabricated in a top-down ap¬ 
proach starting from a (Ga,Mn)As-layer with a Mn- 
content of approximately 5%. The (Ga,Mn)As layer 
we used had a thickness of 15nm and was grown by 
low-temperature molecular beam epitaxy on top of a 
(OOl)-GaAs substrate. In contrast to the experiments of 
Schlapps et alP^ we used as-grown (Ga,Mn)As samples 
without additional annealing before the sample prepara¬ 
tion. First of all, we defined contact-pads for the source- 
and drain-contacts as well as alignment-marks for the 
nanopatterning. This was done using optical lithography 
followed by thermal evaporation of 10 nm Ti and 90 nm 
Au in a standard lift-off technique. After that, the NC 
was defined by means of EBL and subsequent chemically 
assisted ion-beam etching using CI 2 . A two-step EBL- 
process, which allows a precise control of the geometry 
of the nanocontact and a reliable processing, was devel¬ 
oped and is described in Appendix [A] 




FIG. 1: a) Schematic of the PMMA-mask (green / light 
green) defined by a two step EBL-process for etching the NC- 
structure into a (Ga,Mn)As-layer (orange) on top of a semi- 
insulating GaAs-substrate (blue), b) Electron micrograph of 
an NC-device after ion-beam-etching and resist removal. 

The structure of the PMMA mask, used for the two- 
step process, is sketched in Fig. [TJl) . It mainly consists 
of the crosslinked PMMA-line (dark green) of the first, 
high dose (30.000 C/cm) exposure step as well as of a nar¬ 
row gap line from the second, usual exposure step, which 
separates the (Ga,Mn)As layer into two parts used as 
source- and drain-contacts. The two parts are connected 
with each other only at the NC, where the lines of the two 
exposure steps cross each other. This procedure allows 
us to define the width as well as the length of the NC 
by two single lines within independent exposure steps. 


This completely rules out the inter-proximity-effect be¬ 
tween different exposed elements and reduces the mini¬ 
mum size of the NC to the smallest achievable linewidth 
of the two EBL steps. Compared to a single step process 
our approach is robust with respect to minor electron 
dose variations and thus well reproducible. Because of 
this, we were able to fabricate a large number of com¬ 
parable devices and even to control the geometry of the 
NC with a precision of a few nanometers. Fig. 0 >) shows 
an electron micrograph of the central part of a typical 
NC-device taken after the chemically assisted ion beam 
etching and resist removal using a low energy oxygen- 
plasma. After the nanopatterning we covered the whole 
sample with a 30nm thick AI 2 03 -layer grown by a low 
temperature atomic layer deposition process at a tem¬ 
perature of 90° C. The AI 2 03 -layer acts on the one hand 
as the gate-dielectric and on the other hand it protects 
the tiny NC against oxidation. The top-gate contact was 
defined by optical lithography and covers not only the NC 
but also the center part of the whole device. It consists, 
similarly to the source- and drain-contacts, of a 10/90nm 
thick Ti/Au-stack evaporated thermally and structured 
using a standard lift-off technique. 

An effective way to influence the transport behavior 
is to apply an annealing step after the nanopatterning. 
We used an annealing temperature of 150° C and dura¬ 
tions from 30min to 3h. The post patterning annealing 
removes probably some of the defects induced by chemi¬ 
cally assisted ion-beam etching. This can change an ini¬ 
tially insulating sample to one in which Coulomb effects 
prevail or even to a c onducting one. Annealing before the 
nanopatterning 11 ^ 1 , which removes defects induced dur¬ 
ing low temperature molecular beam epitaxy growth, is 
less effective than the post patterning annealing. Hence, 
the intrinsic structure of the NC is dominated by defects 
induced during the nanopatterning rather than by defects 
stemming from the low temperature molecular beam epi¬ 
taxy growth. 


III. MEASUREMENT SETUP 

All low temperature measurements presented in this 
work were carried out at a temperature of about 25mK 
using a 3 He/ 4 He-dilution fridge, equipped with a super¬ 
conducting coil magnet. In combination with a rotatable 
sample holder, we were able to apply magnetic fields up 
to 19T in any direction parallel to the sample plane. In 
order to saturate the magnetization of the device and to 
fix its direction, we applied a constant in-plane magnetic 
field with a magnitude of IT along one of the easy axes of 
the extended (Ga,Mn)As layer. This leads to a situation 
as sketched in Fig. mo- The electrical transport experi¬ 
ments were carried out in a two terminal setup. We per¬ 
formed ac and dc measurements simultaneously by apply¬ 
ing a dc bias-voltage V)j c modulated with a small oscil¬ 
lating ac component V ac . The current I flowing through 
the device was measured using a current amplifier which 
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also converts the current into a corresponding voltage 
signal. The dc measurement using a digital multimeter 
provides the well known I-V^ c characteristic, while the ac 
measurement using a lock-in amplifier offers the differen¬ 
tial conductance G = dI/dV ac of the device. Our device 
could be tuned additionally by an external dc voltage 
(V g ) applied to the top-gate electrode of the device. 


IV. EXPERIMENTAL RESULTS 
A. Room temperature properties 

As mentioned in the introduction, all nanoconstricted 
(Ga,Mn)As-devices investigated in previous studies have 
shown a rath er complex and irregular Coulomb diamond 
pattenf® 1 . This has been explained by assuming that 
several metallic islands are involved in transport across 
the NC. Hence, shrinking the size of the NC should re¬ 
duce the number of islands within the NC and bring 
up a more regular Coulomb diamond pattern. Looking 
for such samples, we investigated many different devices 
with widths and lengths of the NC ranging from lOnm 
to lOOnm. Our experiments revealed that the transport 
properties of these devices are very sensitive to the width 
w of the NC while its length L has only a minor influ¬ 
ence. Wider samples (w > 25nm) show a mainly ohmic 
behavior while the most narrow ones (w < 15nm) are 
fully insulating. Only samples with intermediate widths 
of 15 — 25nm show the typical SET-like behavior, dis¬ 
cussed below. In many cases the room temperature resis¬ 
tance Rnc of the nanocontact already indicates whether 
the constriction is insulating, in the Coulomb blockade 
regime, or ohmic: For R^c/Rs values (with the sheet 
resistance of R s ~ 4kO at 4.2 K) between 10 and 15 
the constriction was in most cases in the Coulomb block¬ 
ade regime for this specific material. However, similar 
to the earlier experiments, all of our SET-like samples, 
even the shortest and narrowest ones, have shown, on a 
first glance, an irregular Coulomb diamond pattern. Be¬ 
low we discuss in more detail transport in the Coulomb 
blockade regime. 


B. Coulomb blockade regime 

In Fig. [2] we present a highly resolved stability dia¬ 
gram of one of our NC-devices in the SET regime. The 
first impression is that the Coulomb diamond pattern is 
very irregular and exhibits frequent vertical discontinu¬ 
ities. Three of them are highlighted by white lines. These 
abrupt shifts can be assigned to charging or discharging 
of local traps in close vicinity to the NC, which, with 
their electrostatic potential, act as local gates. Their ef¬ 
fect can thus be described as an abrupt jump along the 
gate voltage axis. This observation suggests a method 
to reconstruct the stability diagrams with unperturbed 
Coulomb diamonds. We cut the dataset in Fig. [2] along 



FIG. 2: Differential conductance as a function of the bias- 
and gate-voltage of the NC-device in Fig. |T] The measurement 
was done at a temperature of T = 25 mK. A partial irregular 
Coulomb diamond pattern with frequently occurring vertical 
discontinuities is observed. Three of those discontinuities are 
marked by white lines. Cutting the dataset between two of 
these lines gives an undisturbed segment; stitching neighbor¬ 
ing segments together as described in the text and shown in 
the upper inset allows to reconstruct the Coulomb diamond 
spectrum over a larger gate voltage range. 


the white lines and shift the segments on the U g -axis un¬ 
til the diamonds fit onto each other. An example of this 
procedure is shown in the top inset of Fig. [2] In this way 
we obtain, for some parts of the V^-scale, Coulomb di¬ 
amonds which are essentially cleared of potential jumps 
due to charge fluctuations in local traps. The dataset 
displayed in Fig. [3] has been reconstructed from the data 
shown in Fig. [2] and represents the starting point of our 
more detailed analysis. 

The stability diagram shown in Fig. [3] presents char¬ 
acteristic features typical for metallic single electron 
transistors 121 ® but also several anomalies. As expected, 
a series of diamonds of exponentially low differential con¬ 
ductance (black regions with fixed particle number) are 
surrounded by ridges of high conductance. Moreover, by 
further increasing the bias, the differential conductance 
does not drop to zero, see e.g. Fig. §>), allowing to ex¬ 
clude the single particle energy quantization typical for 
quantum dots. Unexpectedly, though, i) the size and the 
shape of the Coulomb diamonds is not regular, ii) some 
of the diamonds are not closing at zero bias (e.g. corners 
between diamond 1 and 2 or between diamond 2 and 3 
as seen from the gate trace in Fig. |3j:)). 

Concerning the first anomaly, it is striking that all 
the diamonds exhibit an individual height as well as an 
individual width. Additionally the diamond labeled 1 
and the diamond labeled 3 are asymmetric: according to 
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FIG. 3: a) Differential conductance of the NC-device of Fig.[2] 
vs. applied gate- and bias-voltages after reconstruction. Dia¬ 
monds labeled 0 to 4 can clearly be identified, b) Differential 
conductance as a function of the bias voltage corresponding 
to the vertical dashed line in a), c) Conductance at 14 = 0 as 
a function of the gate voltage corresponding to the horizontal 
dashed line in a). It shows a conductance peak at the 0-1, and 
a blockade at the other charge degeneracy points, including 
point P. 


the classical orthodox theory®, one would expect that 
all Coulomb diamonds associated to a single island have 
the same size and shape, and that opposing edges of a 
Coulomb diamond were parallel. In the orthodox pic¬ 
ture the two different slopes of a Coulomb diamond are 
related to the capacitive coupling of the island to the 
source- (C s ) and drain-leads (Cu), as well as to the gate 
electrode (Cg). Assuming C g <C C s ,d, the slope of the 
source-line is given by C g /Cd while the slope of the drain¬ 
line is given by —C g /C s (see Fig. 4|. In our case only the 
diamond numbered 2 has parallelsource and drain lines. 
The diamonds labeled 1 and 3, however, exhibit four dif¬ 
ferent slopes, so that we would extract from each two 
different values for C s and Cd or two different values for 
C g , respectively. This suggests that our NC consists ac¬ 
tually of two metallic islands producing a set of nested 




FIG. 4: a) Schematic of a double island structure in a parallel 
configuration. Transport from source to drain is carried by 
two subsequent direct tunneling-processes involving only one 
of the islands. The two islands are characterized by a different 
capacitive coupling to the leads (C B i,Cdi ) as well as by a 
different gate capacitance (C g i) with 1 = 1,2. b) Schematic to 
illustrate the parameter extraction from a regular Coulomb- 
diamond in the framework of the orthodox model, c) The two 
Coulomb-diamonds (ABCD and EFGH) used to extract the 
parameters, are marked by white dotted lines. 


diamonds. 

Fig. [4^t) shows a simple schematic to illustrate our in¬ 
terpretation: the two islands are arranged in parallel, so 
that an electron can tunnel from the source-lead directly 
to each of the two islands and from there in a subse¬ 
quent tunneling process directly to the drain-lead. By 
taking into account the slopes of the diamond edges as 
well as the distance between neighboring charge degener¬ 
acy points we can obtain two different sets of parameters 
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TABLE I: Parameters for the small and large Coulomb dia¬ 
monds (CD) extracted from Fig. §r) assuming a double-island 
structure in the framework of the orthodox theory. The charg¬ 
ing energy U = e 2 /Cs, with Cs = C s + Cd + C g being the 
total capacitance, is also given for reference. 



small CD (ABCD) 

large CD (EFGH) 

Cd 

5.6 x 10" 18 F 

3.0 x 10~ 1S F 

c s 

8.4 x 1CT 18 F 

4.2 x 10" 18 F 

Cg 

28 x 10 _20 F 

9 x 10“ 2o F 

u 

11.2 x 10~ 3 eV 

21.9 x 10“ 3 eV 


(C s , Cu, C g ) from our experimental data. Each set of 
parameters characterizes one of the two islands. One set 
can be extracted from the regularly shaped diamond 2. 
For the other one, we have to reconstruct a second regular 
Coulomb diamond by extending the outer edges of dia¬ 
mond 1 and 3 until they cross each other, see Fig. 4J:). 
The extracted parameters are summarized in Table |T| 
Our analysis is limited to certain gate voltage ranges. 
We attribute this limitation to possible differences in the 
shape and even in the number of participating island as¬ 
sociated to different gate voltage regions. Nevertheless, 
the simple orthodox model gives already a satisfactory 
agreement between experimental and theoretical dl / dV in¬ 
stability diagrams and suggests that transport, in this 
gate voltage range, occurs primarily in parallel across 
two islands of different size in the reconstructed gate volt¬ 
age segments. However, the model presented so far can 
not account for the second anomaly, i.e. a pronounced 
transport blocking observed in the vicinity of the charge 
degeneracy point between the diamonds 1-2 and 2-3, see 
also Fig. §>)• On the other hand, the gap is not present 
at the charge degeneracy point 0-1 and is barely visible 
at 3-4, see also Fig. [3jc). Hence, the gap is assigned to the 
island with the smaller charging energy. In order to ac¬ 
count for this experimental observation, we resort below 
to a more sophisticated transport model that includes 
the ferromagnetic nature of the material. 


V. THEORETICAL MODELING 

In this section we extend the orthodox theory of 
Coulomb blockadd^^ in order to account for the fer¬ 
romagnetic properties of the (Ga,Mn)As samples. Al¬ 
though transport through magnetic islands has been ad¬ 
dressed in the literature) 1 ^ scarce consideration has been 
given, to our knowledge, to the role played by an en¬ 
ergy dependent density of states in the metallic islands. 
The latter, instead, is crucial to explain the anomalous 
current blocking observed in the present experiment. 

To this end we assume that both leads and the metallic 
islands are spin polarized. Fig. [5^,) shows a sketch of the 
magnetization directions expected in the experiments. 
The magnetization of the ferromagnetic (Ga,Mn)As leads 
is rather weak, and can be tuned by an external magnetic 



b) Island 1 Island 2 



M2 


FIG. 5: a) Sketch of the magnetization direction of the leads 
(m s/d ) and of the islands (m/). The magnetization of the 
leads is determined by the direction of the external magnetic 
field. In the constriction, on the other hand, strain effects 
are dominating and the magnetization direction lies parallel 
to the constriction axis. In our experiment the angle between 
the two magnetizations is approximately 6 = |7r. b) Sketch 
of the density of states of the two metallic islands, with the 
spins aligned along the magnetization of the constriction. 


field. It forms in our experiment an angle of 45° (easy 
direction) with the transport direction, set by the longi¬ 
tudinal axis of the NC (z-axis, cf. Fig. [5^)). In the con¬ 
striction, however, the spin polarization axis is strongly 
influenced by strain effects and is expected to be along 
the NC longitudinal axis. 

In order to explain the blockade effects we claim that 
the angle 6 between the leads and the constrictions mag¬ 
netization lies in the range ^ir < 6 < |7r. In other 
words current suppression originates from the fact, that 
the majority spin carriers in the islands and in the leads 
have effectively the opposite polarization. Since only one 
of the two superimposed Coulomb diamond structures 
shows a noteworthy blockade effect, we conclude, within 
our model, that the structure with the blockade stems 
from transport through a fully polarized island, while 
the second island is only partially polarized. 

We describe the islands polarization with an upward 
shift in energy of the minority spin band with respect 
to the majority spin band, see Fig. i>)- The electro¬ 
chemical potential is the external parameter which deter¬ 
mines whether the island is partially or fully polarized. 
Partial polarization is obtained if the chemical potential 
Ha (a = 1, 2) lies above the bottom of the minority spin 
band, full polarization when the chemical potential lies 
between the bottom of the majority and of the minority 
spin bands. 
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In our model the tunneling of a source electron of the 
majority spin species (conventionally the spin up) to a 
fully down polarized island is highly suppressed for low 
bias voltages since no spin up states are available near 
the Fermi level. For bias voltages which are large enough 
to access also the minority spin band (age 14 > B+, cf. 
Fig- ID 3 )), the suppression is lifted and an increase of the 
current is expected. For the partially polarized island 
both spin species can be accessed already at the Fermi 
energy and no suppression is observed. 

A. Model Hamiltonian 

We describe the nanoconstriction with a system-bath 
model aimed at mimicking the structure of the two is¬ 
lands contacted to source and drain leads sketched in 
Fig. |4jr). The total Hamiltonian is 

H = H S + H T + H h , ( 1 ) 

where 

Hh = E v ka c\ k <j C v ka (2) 

??e{S,D} ka 

denotes the Hamiltonian of the two spin polarized leads. 
We assume to have a flat, but spin dependent, density of 
states (cr =t / |) 

At = Art = ^ A, (3) 

which depends on the polarization p v of the leads (—1 < 
Pri < 1)- The metallic islands (a £ {1,2}) in the 
nanoconstriction are modeled by 

Hq — ^ i ^ ) Cckzt d a j T d ocir Cl.. ( \ 1 -V,, 

aS{l,2} it m 

+ U ^N a (N a - 1 )|, 

and have in general a different spin quantization axis as 
the contacts. We define r = ±1 for spin +/—, respec¬ 
tively, using the spin-quantization axis of the nanocon¬ 
striction. As already mentioned, we account for the ferro¬ 
magnetic properties of the metallic islands by assigning 
spin dependent energy levels, e a i T , and consequently a 
relative shift of the density of states for the two spin di¬ 
rections, A ex (Fig. [5^)). The long range Coulomb inter¬ 
actions are included within a constant interaction model, 
where U a is the charging energy of the island a. The 
effective coupling of the gate electrode to the metallic is¬ 
lands is taken into account by the term proportional to 
a g eV g , with a g = C g /Cs being an effective gate coupling 
parameter and V g the gate voltage. The two metallic is¬ 
lands and the leads are weakly coupled by the tunneling 
Hamiltonian 

= EE ( J'TjOicr'U'aT (®) ^rikcr daiT + h.C.'j, (g) 

ioLT r]k(T 


where we defined the function Uf+(0) = u±-(6) = 
cos(0/2), Uf-(0) = u±+(9) = *sin(0/2). It results from 
the non-collinear spin quantization axis of the islands and 
the leads. Since the two axis are rotated by an angle of 
9 in the y-z-plane with respect to each other, the trans¬ 
formation conserves the spin during tunneling. 

B. Density of states of the metallic islands 

Some of the experimental observations can only be un¬ 
derstood if the energy dependence of the density of states, 
in particular the presence of different band edges for mi¬ 
nority and majority spins, is accounted for. Specifically, 
we define the spin-dependent density of states of island 
a as: 

ffcr(e) = 

= den Q(e + W — t A e x/ 2 ) 0 (VF -I- rA e x/2 — e) ( 6 ) 
~ 9ut f (e + W — rAex/2), 

where W is the spin independent contribution to the 
bandwidth, and A e x the exchange band splitting of the 
ferromagnetic metallic island. The parameter g aT de¬ 
fines the strength of the density of sates. Since the W is 
the largest energy scale considered in the following, the 
upper limit of the density of states can be set to infin¬ 
ity. In the last line of Eq. © we have approximated the 
left Heaviside function by f~ = 1 — f + , with f + the 
Fermi function; this allows us to further proceed analyt¬ 
ically in the calculation of the transport properties. The 
density of states is also sketched for clarity in Fig. ©>)• 
For later reference we define 13“ as the energy difference 
between the bottom of the band of the corresponding 
spin species r and the chemical potential of the island a: 
= -W + rAex/2 - \x a - 

C. Transport theory 

In the following we briefly outline the main steps 
leading to the evaluation of the transport characteris¬ 
tics, emphasizing the new ingredients entering our trans¬ 
port theory. For more details we refer to the Appendix 
|B| The framework is the orthodox theory of Coulomb 
blockade 1 ®^, extended to the case of ferromagnetic 
contacts 17 ^ and valid also for fully spin polarized metallic 
islands. The explicit derivation of the tunnelling rates 
should illustrate the crucial role played in our theory by 
the energy dependent density of states. 

The theory is based on a master equation for the re¬ 
duced density matrix of the islands, up to second order 
in the tunneling Hamiltonian. Since the two metallic is¬ 
lands are assumed not to interact with each other, the 
corresponding density matrices obey independent equa¬ 
tions of motion (see Appendix [B]). Moreover, the metallic 
islands are assumed large enough to posses a quasi con¬ 
tinuous single-particle spectrum, but small enough that 
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their charging energy dominates the tunnelling processes 
that change their particle number. We further assume 
that, in between two tunnelling events, the islands relax 
to a local thermal equilibrium. Under these assumptions 
the reduced density matrix of island a can be written as 

„ f e -/3fls,c 't 

PredW = E 1 ' Pn <* “W- f P N a (t), (7) 

N a ^ *N a J 

where Hs a is the part of the system Hamiltonian asso¬ 
ciated to the island a, Vn„ is the projection operator on 
the N a -particle subspace and Z^ a = Trs(7 ' > N a e~ l3Hs ’ a ) 
is the corresponding (canonical) partition function. By 
projecting the master equation on the iV a -particle sub¬ 
space and tracing over the islands degrees of freedom, we 
keep only the occupation probabilities Pn„ of finding the 
island occupied by N a electrons as dynamical variables. 
In the stationary limit we find (see Appendix |b|) 


Trs {V N j° 0 o} = 0 



Eventually, the stationary current through lead ?y reads 

u = - e E E { K^ Na+1 - W ■ (9) 

olo N a ^ ' 

In Eqs. ([Sf and § the rates are defined as 
r ^ +1 ^“ = E Kr(<?)| 2 b~(AE% a - a v eV b ) 

T 

x |e(AU^ + - a^eVb) - F(B ?)|, 

r$? N - +1 = E \^r{9)\ 2 b+(AE% a - a v eV h ) 

T 

x |f(B“) - F(AE% a + B° - c^eW) j, 

(10) 

and are expressed in terms of the normal state resis¬ 
tance = h/{2 / ne 2 \t r}a<T \ 2 g aT D r) ) and the functions 

b ± (x) = l/(e ±/3x — 1) and F(x) = x/(eP x — 1), with 
/J = \/{k^T) the inverse temperature. We account for 

the asymmetric bias drop with the bias coupling con- 

C ~\~C /2 

stants defined as «s/d = ± d/ c s S —• Further, we de¬ 
fined the grand canonical addition energy 

A £jv Q = a g e Vg + U a N a + / i a — fj, o 

= (pN a + l — MO (N a + 1)) — (E/v q — MO Not) 

which must be paid in order to increase the electron num¬ 
ber on island a from N a —► N a + 1. We denote Mo the 
chemical potential of the leads at bias 14 = 0. 


The rates given in Eq. ( [To] ) differ from the ones of 
the orthodox theory of Coulomb Blockad^ 12 ^ even in 
their spin dependent variation 1 ' due to the energy de¬ 
pendent density of states and the explicit dependence on 
the band edges. The latter introduce a new source of 
current suppression associated to the absence of states 
with a specific spin species. These rates represent the 
main theoretical contribution of the present work. For 
the chemical potential lying far above the bottom of the 
bands, the theory recovers again the limit of the classical 
orthodox theory of Coulomb blockade. Namely, in the 
limit B —> —oo: 

lim ±b ± (x){F(B)~ F(x + B)\ = F(±x). (12) 

B—> — oo l ’ 


VI. THEORETICAL RESULTS 
A. Comparison with the experiments 

The results of our simulation are reported in Fig. [6^,), 
with the differential conductance shown as a function of 
the bias and gate voltage. We see the same nested dia¬ 
mond structure as in the experiments. In our theory the 
diamonds at the charge degeneracy points labeled 0-1 and 
3-4 close. Between the diamonds 1-2 and 2-3 the differ¬ 
ential conductance is suppressed for bias voltages smaller 
than a certain threshold bias. Fig. [6jo) shows a bias trace 
calculated at the charge degeneracy point 1-2, for two 
different angles 6 between the magnetization vectors of 
the leads, m a , and the metallic islands, m/. It shows a 
suppression of the differential conductance at point (P) 
with respect to point (Q). The width of the suppression 
region corresponds to the one observed experimentally in 
Fig. [3 Jd) and is proportional to B +, the energy difference 
between the bottom of the minority spin band and the 
chemical potential of island 1 (cf. Fig. [5 ]d)). In contrast 
to the experiments no full blockade can be observed at 
(P). A change of the orientation of the magnetization di¬ 
rections from 0 = 7r (dashed red line) to 9 = j7r (solid 
blue line) is shifting the curve upwards. Besides the con¬ 
stant shift the two curves are qualitatively the same. 

To emphasize the effect of the islands degree of po¬ 
larization on the suppression mechanism, a conductance 
trace at 14 = 0 of a full polarized island 1 is compared 
to the case of a partial polarized island 1 in Fig. I/O- 
Partial polarization is achieved by shifting the electro¬ 
chemical potential of island 1 by 12 meV up in energy. 
The solid blue line shows the full polarized case, where 
the two larger peaks correspond to the larger Coulomb 
diamond (island 2). The peak observed in the experi¬ 
ment (Fig. [3j:)) we ascribe to transport across this par¬ 
tially polarized island. Although the theoretically pre¬ 
dicted second peak is missing in Fig. [3j:) we note that 
the corresponding blockade between diamond 3 and 4 is 
much less pronounced than between, e.g. 2 and 3. This 
asymmetry between the degeneracy points 0-1 and 3-4, 
however, cannot be accounted for by our model which 
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FIG. 6: a) Calculated differential conductance of two spin po¬ 
larized metallic islands with spin polarized leads. The island 
with the larger charging energy is assumed to be partially po¬ 
larized, while the island with the smaller charging energy is 
fully polarized, b) Bias trace through a charge degeneracy 
point of the fully polarized island. The gate position of the 
line trace is marked as a dashed line in a), c) Gate traces at 
Vb = 0 for a full (solid blue line) and a partial polarization 
(dashed grey line) of island 1. Island 2 remains partially po¬ 
larized. In the fully polarized cased the conductance peaks of 
island 1 are suppressed with respect to the partially polarized 
case. For island 2 both curves are identical. The parame¬ 
ters used to obtain this figure are: «si = 0.4, 032 = 0.42, 
U\ = 11.2 meV and U 2 = 21.9 meV in accordance with the 
parameters for the capacitive couplings of Table [T] Moreover 
RT t = 0.57-10 3 h/e 2 , = 1.4-10 3 /i/e 2 , B\ = 2meV, B x _ = 

— 18meV, B\ = — lOmeV, B 2 = — 35meV, pi = —42meV, 
fj ,2 = —32meV, p v = 0.8, and k^T = 0.07 meV. For the full 
polarized island 1 in c) B\ = — lOmeV, Bl = — 30meV. 


predicts a periodicity of the Coulomb oscillation pattern. 
The four smaller peaks in Fig. [6)c) belong to the smaller 
Coulomb diamond structure, corresponding to island 1, 
i.e. the fully polarized one (Fig. |5 )d)). Even though the 
conductance is not completely suppressed as in the ex- 
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FIG. 7: All tunnelling rates for the island 1 plotted at a 
charge degeneracy point as a function of the bias voltage Vb- 
The angle between the two magnetization direction is 9 = tv. 


periment, the conductance peaks are strongly reduced 
with respect to the partially polarized case. In the lat¬ 
ter (dashed grey lines) no suppression is present and the 
conductance peaks of island 1 are by a factor of 4 larger. 
Below we address a possible reason for the incomplete 
blocking within the model. Since the parameters of is¬ 
land 2 are kept the same, both for the fully and partially 
polarized cases the corresponding conductance peaks are 
not changing. Despite the fact that a comparison of the 
calculated gate trace to the experimental one in Fig. [ 3 J 2 ) 
reveals some limitations of the model, the essential fea¬ 
ture, i.e. the suppression inside the large Coulomb dia¬ 
mond, is reproduced. 


B. The mechanism of current suppression 


For a better understanding of the mechanism under¬ 
lying the blockade, we derive analytically the differential 
conductance for the island 1 at the two points (P) and 
(Q) marked in Fig. Do). For simplicity the case 8 = tt is 
considered, since qualitatively the blockade mechanism is 
the same in both cases. 

Notice that both P and Q correspond to a gate volt¬ 
age such that A = 0, i.e. at the charge degeneracy 
point of the N-N + 1 transition. To obtain the differen¬ 
tial conductance, according to Eq. ([8]) and ([9]), the tran¬ 
sition rates are required. For simplicity we have 

dropped the subscript 1 from the excitation energy A 
since we will refer from now on always to the same island. 

In Fig. [7] we show the transition rates as a function 
of the bias 14- To simplify the notation, we replaced 
r vc^ N±1 Notice their linear dependence on the 

bias above a certain threshold. Thus, in that bias range 
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one can approximate them as: 


I&j. = -%Vb, 

r DJ. = 

= ^Dt + ®Dt 14 j 


(13) 


where ^Df is a constant accounting for the threshold 
bias, and 


_ 27 re i + 075 2 

^770- — j -2 -^0 2 Qcr^T] ^77 1 • 


(14) 


Here, Ho = is assumed to be independent of the 
lead. For the point (P) within the first plateau only 
the rates with a =J,, namely Fj^ and Tg^, are nonzero. 
Hence, according to the principle of detailed balance: 
T Dl p N = F<Hv +1 . Imposing probability conservation 
we find P/v = Fg|/(r^ + FgjJ. Thus the stationary cur¬ 
rent equals 4 P) = -er^P/v = -eT> ; r^/(T^+r^) <x 
(1 — p) 2 , which is suppressed for a large spin polarization 
p. Here the polarization p is assumed to be equal for 
both leads. 

At the point (Q), only one additional rate, rj4, is con¬ 
tributing (the rate r)4 is zero due to the lower bound 
of the density of states). In this bias range the equa¬ 
tions of detailed balance and probability conservation 
yield Py = r|^/(r^ + + F^). The resulting sta¬ 

tionary current is then 1^ = — e(Fj^ + F^)P.y = 

—e (^Dj. + + ^Df + r<) oc (1 - p). Again 

the current is suppressed for large spin polarization. 

Inserting Eq. (131 into the current expressions at the 


points P and Q we find 


t(p) _ B sil3ni T/ 

*D ~ e K> n *4 

^ Hr 


->S 4 . 


(15) 


and 


j(Q) = c +BsiAutVb - Hs;(Hdt + BpiWb ^ 


Ad| + (#Dt + - Ssj.)Vb 


Taking the ratio of the two differential conductance 
plateaus, i.e. the ratio of Eqs. (Dl) and (D2), we find 


dldl^ 

a± D jUi D 


1 


dV -b 1 dV b odI^dI 2 — asl^sl 2 

2 _ U_.|2 (i-p)h 


x anltuf - as|ts|' 


(l+p)p t + (l-p)p| 


(17) 


Thus, within our simple model, the ratio TZ of the height 
of the two plateaus is limited by 


omltsl" 


odI^dI 2 - cks l*sl 


<n< 1. 


2 — — 


(18) 


In other words the ratio of the two differential conduc¬ 
tance plateaus is limited in our theory, leading to some 
discrepancy with the experimentally observed ratio, cf. 


points (P) and (Q) marked in Fig. [3 ]d). Since the param¬ 
eters a v are determined experimentally, the only possibil¬ 
ity to change the ratio is to modify the coupling constants 
\t v \. However, the increase of the coupling constants nec¬ 
essary to fit the experimental value, would lead to a huge 
asymmetry in the stability diagram which is not observed 
experimentally. Despite the discrepancy between TZ and 
the experimental ratio, we think that the theory clearly 
suggests a mechanism which can lead to a suppression of 
the conductance due to spin polarization in the frame¬ 
work of an orthodox theory of Coulomb blockade. To 
better fit the experiments a more realistic energy depen¬ 
dence of the density of states which also accounts for 
valence bands is necessary. With such an energy depen¬ 
dence the rates can change their slope as a function of 
the bias voltage, leading to an even more pronounced bias 
dependent suppression of the differential conductance. 


VII. CONCLUSION 

In this work we have reported on a detailed study 
of the transport characteristics of nanofabricated nar¬ 
row constrictions in (Ga,Mn)As thin films. By means 
of a two step electron beam lithography technique we 
have fabricated well defined nanoconstrictions of differ¬ 
ent sizes. Depending on channel width and length, for 
a specific material, different low-temperature transport 
regimes have been identified, namely the ohmic regime, 
the single electron tunnelling regime (SET) and a com¬ 
pletely insulating regime. In the SET, complex stability 
diagrams with nested Coulomb diamonds and anomalous 
conductance suppression in the vicinity of charge degen¬ 
eracy points have been measured. 

In order to rationalize these observations we proposed, 
for a specific nanoconstriction, a model consisting of two 
ferromagnetic islands coupled to ferromagnetic leads. In 
particular, the angle 8 between the leads and the islands 
magnetization lies in the range < 8 < r. Moreover, 
the full polarization of one of the metallic islands is cru¬ 
cial. We studied the transport characteristics of the sys¬ 
tem in terms of a modified orthodox theory of Coulomb 
blockade which takes into account the energy dependence 
of the density of states in the metallic islands. The latter 
represents an important generalization of existing formu¬ 
lations and is determinant for the qualitative understand¬ 
ing of the present experiments. In fact, the explicit ap¬ 
pearance of the minority spin band edge in the expression 
of the tunnelling rates yields a pronounced conductance 
suppression at the charge degeneracy points. To account 
for the full suppression of conductance observed in the 
experiments the simple model used in this work should 
be further improved. For example the hole character of 
the charge carriers and associated spin orbit coupling ef¬ 
fects are not captured by our model. Furthermore, it is 
straightforward to combine the present theory with mi¬ 
croscopic models that allow for a realistic description of 
the islands density of states. 
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Appendix A: Sample fabrication 


Since the leads are considered thermal baths of noninter¬ 
acting fermions, /5 l reads 

p —m 

Pl =-=-■ (B3) 

^L.G 


1. Two-step EBL fabrication process 

Both steps are based on the standard EBL resist 
poly-methyl-methacrylate (PMMA). In the first step one 
exposes the resist using an extremely high line-dose 
(approx. 30.000pC/cm) in order to define a narrow 
crosslinked PMMA-line. This line is very robust and does 
not get removed by common organic solvents like acetone. 
Hence, after cleaning the sample in a bath of acetone, 
the crosslinked PMMA-line remains on top of the sample 
while the unexposed PMMA is removed from the sample 
surface. For the second step the sample is again coated 
with a fresh layer of PMMA resist. This time one uses 
a common dose (approx. 2000pC/cm) in order to ex¬ 
pose a second line perpendicular to the crosslinked one. 
After removing the exposed resist using a standard devel¬ 
oper solution consisting of isopropyl alcohol and Methyl- 
isobutyl-ketone (MIBK), we get the patterned mask for 
the subsequent ion-beam-etching, shown in Fig. |iJl) . 


Appendix B: Equation of motion for a orthodox 
theory of Coulomb blockade 


Further, we assume that due to fast relaxation processes 
in the leads, the density matrix can be written as pi(t) = 
/5 re d,i(i)/5L + 0(H T ) , with /5 red ,i = Tr Ly o. Moreover, due 
to the independence of the two metallic islands p re d(t) = 
/5(ed(^)pred(t) and each component obeys the following 
equation of motion: 




= -^[H S ,p? ed (t)] 


1 

¥ 



Pr ed (t)p L \ 


(B4) 


where a — 1,2 labels the metallic island. 

For the system we assume that the metallic islands are 
large enough to posses a quasi continuous single-particle 
spectrum, but small enough that their charging energy 
dominates the tunnelling processes that change their par¬ 
ticle number. Furthermore, it is assumed that the islands 
will relax to a local thermal equilibrium on a time scale 
shorter than the inverse of the average electronic tun¬ 
nelling rate. Under these assumptions, the reduced den¬ 
sity matrix can be written as 


In this appendix we derive an extension of the orthodox 
theory of Coulomb blockade for the case of spin polarized 
contacts as well as of a spin polarized metallic island. In 
particular we will consider explicitly the lower bound of 
the density of states in the metallic island. 

The transport theory is based on the Liouville-von 
Neumann equation for the reduced density matrix in the 
interaction picture 

(Bl) 

which we expand to second order in the tunneling Hamil¬ 
tonian Ht■ Prior to t = 0 the system and the leads do 
not interact and the density matrix can be written as a 
tensor product of the density matrices of the subsystems 

(B2) 


Pr“edW = E^ 


N a 


o-pHs,, 


■>N a 


P N a (t), 


(B5) 


with Z/v q = Trg {'Py ct e and 

P N a — ^ ^ ; 

{rule, 

ni-N a 


(B6) 


is the projection operator on the _/V a -particle subspace. 


Notice that in Eq. (B5), due to the projector opera¬ 


tor Vn 0 > the only statistically relevant term of the sys¬ 
tem Hamiltonian H s , a is hg = Yia £ aiadl ia d aia . The 
term e -P(^ LN tx(N 0 ,-i)+a e .eV s N a ) b ecomes a constant and 
is canceling out in the density matrix. Inserting explicitly 
in Eq. (B41, we find 


P = Ps( 0) ® p L = /3 s(0)/3l 
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Trs {7Vj“dW} 


- 'y^. 'y^. trja<r'U' ( j T (d) t^i aa iU a i T , (0) J dt / 

r)i)' kiar ® ^ 

iii / / 

k i cr r 

Trs | V Na dLr ^^./(-t'OPredWj Tl'L { C vka 
T Trs | V Na dair dL/ T /,i(-^")ft“ed(i)| Tl 'L { e; te C^fe'a'.K-i'OPL} 

- Trs | dai'r',1 {-t")V Na dLr PredW} Tr L { c I'k'c'ji-t") Cvk* PL } 

- Trs | dL' T ',I(-^")^JV Q <W Pre d w| Tr L { VfcV'.lM") cj, fe(T PL } 


+C.C. >. 


(B7) 


r 


In the following we are analyzing the first term of Eq.(B7) 


in more detail, the other terms can be evaluated in com¬ 
plete analogy. The calculation of the trace over the lead 
degrees of freedom gives 


freedom another approximation is necessary. By taking 
the average in the grand canonical ensemble, the particle 
number is determined by the chemical potential and we 
can remove the projection operator: 


Tl'L { C^fccr C^/ fc / £r / I ( t )/5l} 


= g * ) f 


f (Pr)k Pr])^kk' & r/rj 1 $(7 Ct' i 


(B8) 


where the time evolution of the creation and annihilation 
operators of the leads is given by c' vkal (t) = e* E ’ lt c' ka - 
For the system operators the time evolution can be car¬ 
ried out in a similar way, keeping in mind that the parts 
proportional to the total number operator can be factor¬ 
ized 


Tr s \ V Na d aia d aiV ,i(-f")p r “dW 


= Trs < VN a d aia d a iv',i(—O 


'P: 


J N a 


N a 


f „ + - e~P\ n S 

Trs \ d aiCT d Qi / CT Li 






Ma, N a 


Pi 


N a 


(BIO) 


Tr s <i PN a dl iT d ai / r qi (~t")p? ed (t) 


X Trs i PN a d aiT d ai'r‘ 


P?ed(t) y 


(B9) 


In order to perform the trace over the system degrees of 


This approximation becomes exact in the limit of 
N oo. In presence of a quasi-continuous energy spec¬ 
trum of the islands we can further drop the N a depen¬ 
dence of the chemical potential, for small relative varia¬ 
tions of N„. 


The trace in Eq. (B9) can now be evaluated in the 


standard way and it yields Fermi functions. Inserting 
the results for the traces in Eq. (B7) we obtain: 


Tr s 


{P»AWl = -pEI i* 


l 1)OL(T I 


-( 0 ) 1 2 


dt h 


i) kicrr 


e i(-E vk+ e^ +as eV s +U a (N a -!))*"/+( Cq . t _ p a )f-{E„ k - p v ) P Na (t) 

+e -U~^ k +e air + as eV g +UN a )t" f -( Ca . T _ /Xa )/+( jE??fc - P N Jt) 


_ E vk + 


r +a s eV s +U(N a -l))t" j- ^ 


fJ"l)') PN a — 1 (0 

_ e -i(-^ k +e ai r+ as eV e +UN a )t" f+ ^ aiT _ - /Z„) P Na + 1 (t) 


+C.C. 


(BH) 
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Since we are only interested in the stationary solution 
of the master equation, we send t —> oo and use the Dirac 
identity 

dte lult = tt5(uj) + * lim Im ( --— j (B12) 

j)->o \uj + ir) J 



to evaluate the integrals. Due to statistical averages no 
coherences are possible in the master equation and the 
two complex conjugated parts can be summed up. We 
find 


Tr s {V N J^} = IVJ 2 KtWI 2 j 

r] kiar ^ 

-\- € a i r + OLgCVg + Ua(Na> 1)) Mck)/ (-^ r/k l^rj) -P/Vq. 

+5( — + e a i T + a g eV^ + UN a ) / ( e a i T — /i a ) f + (E^k — ^rj) PN a (B13) 

6( E^k + C a i T ~\~ OtgCVg ~\~ U(N a 1)) f {tair Ma)/ P]y a — 1 

(i ( E^k p Cazr T CfcgeVjr + UNa) f^~(c a i T (E^k [i T j ) -FjV^ + 1 r ■ 


Further, we consider the continuum limit of the states in 
the quantum dot 

/ OO 

de 3ar(e), (B14) 

z -°° 

with g aT (e) being the energy dependent density of states 
in island a with the spin r, defined in Eq. ©• For the 


leads 

/ OO 

dE D vai (B15) 

K 

where D v „ is the density of states of lead 77 which is 
considered in the flat band limit. The integration over 
the lead degrees of freedom gives: 


Trs {pN a P^,} = 0 = — ^ ^2 IWl 2 \u aT {0)\ 2 D va fdeg aT {e){ 

rfCTT ^ 


f + (e-n a )f (e +AE Na _i - fi^) P Na 
+ / ( e — d-a)f + {(- + ^EN a — f-lrj) Pn a 

-/“(e - H a )f + {e + AE Na -i - n v ) Pn^-i 
-/ + (e - y- a )f~(e + AE Noi -^ v )P Na+ 1|. 


(B16) 


where AEjy a = UN a + a g elg. In a last step we insert 
g ar (e), see Eq. Q in the main text, and the remaining 
integral can be done by using the following identities: 

f + {x)f~{y) = b + (x - y)(f + (y) - f + {x )), (B17) 


dx ( f + {x ) - f + (x + uj)) = uj, 


(B18) 


f dxf + (x + a)f (x + b)f [x + c) = f dx b + (a — b) (f + (x + b) — f + (x + a)j / (x + c) = 
=b + (a — b)l\b + (b — c) J dx ^/ + (a; + c) — f + (x + 6)^ — b + (a — c) J dx ^f + (x + c) — f + (x + a) 
=b+(a - b) (^F{b - c) - F(a - c)^ 


(B19) 
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b ± (x ) a nd F(x) are defined in the main text just below 
Eq. (flO]) . Using these identities yields the final result 


Trs {V N J° 0 o} = 0 

E l T'Na—tNa — 1 D ~rN a ^N a +l P 

j r/aa 9 N a ^ rjacr 

r)cr ^ 

_l_ Y^Na—l—tN P . -pNa + l^-N p I 

+ i / ? Q !0 - + 1 ncxa ^N a +1 >■ 

(B20) 


In the interaction picture the total particle number op¬ 
erator of lead r] 7 N v , is not evolving in time since it 
commutes with the unperturbed part of the Hamiltonian. 
Therefore, the current reads 


Appendix C: Current 

Finally we briefly outline the derivation of the current 
formula. The current is defined as 

U = ej t (N v ) ( t ). (Cl) 


Iri = eTr s+ L \ N v j t Pi(t) 


= - ^ Tr s+ L N v [H T ,i(t), pi( 


~Jp dt ' Tr s+L j N 


H T .i(t), [HTj(t'), pi(t')\ 


(C2) 


where we expand 4_pi(t) up to second order in H + . The we replace P\{t r ) —> pi{t). Exploiting further the cyclic 

first term of Eq. ( |C2| ) vanishes since only a odd number invariance of the trace we find 

of operators appear in the trace. In the second term 


Iv ~ v L dt ’ Tr 


[Nr,, H T ,l{t)] , HT,l(t') 


Pl{t) 


2e R 

= -v** 


dt' Trs+L j [N v , HTj(t)]HT,i(t')pi(t) 


(C3) 


In the last step we exploited the anti-hermiticity of 
[N,,, . Following the same steps as in the deriva¬ 

tion of the master equation, one can identify the rates, 
and one finds the well known expression of the current 


Appendix D: Calculation of the differential 
conductance 


Differentiating Eq. (151 with respect to 14 and insert¬ 


ing the definition of Eq. (14) yields the differential con¬ 


ductance of the first plateau: 


4 = 


e EE 

acr N a 


N a —y N a -f- 1 p _ N a —y N Q — 1 p 


N n 


(C4) 


di£_ 

d(V h ) 


27re 2 _ (1 -p) a D |to| 2 Q:s|4| 2 

h o9i 2 ap|f D | 2 -asl4l 2 ' 


(Dl) 


To calculate the differential conductance at this point we 


differentiate Eq. (16) with respect the bias voltage and 






















14 


find 


dl^ _ 2 ot + 2^7K+W 

dV b 6 ( 7 + W 


(D2) 


where we defined a = /3 = -B$i(Bui + S D |), 

7 = *4Dt> an d d = —Bsi + Bui + Bui- 1 11 order to find 
the value of the differential conductance plateau we have 
to consider the high bias limit and we find 


Inserting back the physical constants we find 


lim dI ° 


(Q) 


2 2 tt (1 -p) . |2 

vwoo dV h =e Y Do9i 2 as|fs| 

x _ ap\tu\ 2 {(l +p)9t + (1 -P)h) _ 

X ~ P)9i a s\ts\ 2 + ((1 +P)9t + (1 - p)9i)oiD\tu\ 2 ' 

(D4) 


lim 

Vb-^-oo 


dI s Q) 

~dK 


2 @_ _ 2 ^St(^Dt+gp;) 

e 5 _e -Bsi + Bui + Bui 


(D3) 
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